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ABSTRACT 

The progenitors of Type la and some core collapse supernovae are thought to be stars in binary 
systems, but little observational evidence exists to confirm the hypothesis. We suggest that the colli- 
sion of the supernova ejecta with its companion star should produce detectable emission in the hours 
and days following the explosion. The interaction occurs at distances ~ 10 11 — 10 13 cm and shocks 
the impacting supernova debris, dissipating kinetic energy and re-heating the gas. Initially, some 
radiation may escape promptly through the evacuated region of the shadowcone, producing a bright 
X-ray (0.1-2 keV) burst lasting minutes to hours with luminosity L ~ 10 44 ergs s _1 . Continuing 
radiative diffusion from deeper layers of shock heated ejecta produces a longer lasting optical/UV 
emission which exceeds the radioactively powered luminosity of the supernova for the first few days 
after the explosion. These signatures are prominent for viewing angles looking down upon the shocked 
region, or about 10% of the time. The properties of the emission provide a straightforward measure of 
the separation distance between the stars and hence (assuming Roche lobe overflow) the companion's 
radius. Current optical and UV data sets likely already constrain red giant companions. By system- 
atically acquiring early time data for many supernovae, it should eventually be possible to empirically 
determine how the parameters of the progenitor system influence the outcome of the explosion. 

Subject headings: 



1. INTRODUCTION 

Observations of supernova light curves and spectra 
have allowed us to characterize the outcome of the ex- 
plosion - the burned and ejected stellar debris - in re- 
markable detail. But we still know very little about the 
starting point. In the most widely considered scenario, 
Type la supernovae (SNe la) result from carbon/oxygen 
white dwarfs which reach a critical mass by accreting 
material from a non-degenerate companion star. Obser- 
vational confirmation of the binary nature of the pro- 
genitor system is lacking, however, and almost nothing 
is known about the properties and diversity of the com- 
panion stars. The proposed progenitors should be rather 
dim, so it is not surprising that we have so far failed 
to fi nd them in pre-explosion images of the host galaxies 
(e.g., Maoz & Mannucci 2008 ). One therefore seeks other 
means of inferring the presence of a stellar companion. 

A few minutes to hours after the supernova eruption, 
the debris ejected in the explosion is expected to over- 
run the companion. The star is shocked by the impact, 
and it's envelope pa rtially stripped and ablated, but it 
survives t he ordeal ([Wheeler et al.||1975| [Fryxell fe Ar- 
nett|1981| |Chugai|1986| |Livne et al.|1992| [Marietta et al. 
2000[ jPakmor et al.||2008j r Observations of the rem- 
nant of Tycho's 1572 supernova turned up a high- v elocity 
G star, claimed to be the runaway companion (Ruiz- 
Lapuente e t al. 200 41). T his identification is still debated 
( |Kerzendorf et al. 20091. Meanwhile, several attempts 
to look for evidence of stripped hy drogen in the spec- 



tra of SNe la h ave detected nothing ( Mattila et al.|2005 
Leonard||2007l. The supernova ejecta is distorted by the 



collision, whic h should lead t o polarization of the super- 
nova light ([Ka sen et al. 12004] ) . But alt hough polarization 
has been detected in several SNe la ( |Wang fc Wheeler] 



2008 1 , there is nothing to unambiguously indicate that 
this asymmetry relates to companion interaction. 

While the previous investigations have focused on the 
long term consequences, one might wonder: could we 
see the collision itself? Here we can draw an interest- 
ing parallel with core collapse supernova explosions, in 
which a shock wave propagates through the envelope of 
a massive star. When that shock front nears the stellar 
surface at radius R « 10 11 — 10 13 cm, the post shock en- 
ergy v ents in an X-ray breakout burst lasting minutes t o 



hours flKlein fc Chevarier|1978~l|Matzner fc McKee|1999 l 
In the days that follow, optical/ UV radiation continues 
to diffuse from the deeper layers of shock heated ejecta. 
Eventually, the energy deposition from radioactive 56 Ni 
decay takes over. The early shock luminosity phase, how- 
ever, has been obser ved in several e vents, e.g., SN 1987A 



flArnett et al.]|l989l and SN 1993J ( | Wheeler et ^[1993 1 
while the shock breakout itself was caught for SN 2008D 
( [Soderberg et al.|2008| |Modjaz et aT1|2009[ ). 
The same physics applies to Type la sup 
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Lype la supernovae, but 
because the radius of the progenitor white dwarf is so 
small (i? w d ~ 2 x 10 8 cm) the breakout emission should 
be extremely brief and the early luminosity remarkably 
dim. The problem is that when energy input occurs 
at small radii, adiabatic losses in the rapidly expand- 
ing (v ss 10 9 cm s _1 ) ejecta are overwhelming, and the 
thermal shock deposited energy is converted to kinetic 
energy on the expansion timescale (i? w( j/w ~ 0.1 sec), 
which is much shorter than the diffusion timescale. 

As it turns out, the separation distance between the 
white dwarf and its companion star is presumed to be a ~ 
10 11 — 10 13 cm, comparable to the radii of core collapse 
progenitors. When the supernova ejecta collides with its 
companion, the impacting layers are re-shocked and the 
kinetic energy partially dissipated. If the geometry is 
favorable, some of this energy might escape straightaway 
in a prompt burst, which will be followed by a longer- 
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Fig. 1. — Hydrodynamical calculation of a Type la supernova 



colliding with a red giant star (R+ 



10 1 



■■ 2.5 X 10 13 cm) 



a: Density structure during the prompt emission phase (t = ti/2). 
The companion star (drawn in blue) diverts the flow and carves a 
hole in the ejecta. The black contour shows the region where the 
shocked ejecta temperature exceeds 3 X 10 5 K. b: Density structure 
at a later phase (t = 2tj, note change of scale). The shell of shocked 
ejecta has expanded to partially refill the hole. The black contour 
shows the region where the temperature exceeds 10 5 K. 



lasting tail of diffusive emission - the analogues of shock 
breakout and its aftermath in core collapse events. In this 
case the emission provides a measure not of the stellar 
radius i? w d, but of the separation distance a. 

Here we develop an analytic description of the colli- 
sion dynamics and subsequent radiation transport which 
suggests that early time observations of supernovae at 
X-ray through optical wavelengths should offer a power- 
ful means of confirming the presence of a companion star 
and constraining its parameters. 

2. COLLISION DYNAMICS 

In the single degenerate scenario of SNe la, the com- 
panion stars are thought to be eit her slightly evolved 



main sequence stars or red giants (Branch et al 
Hachisu et al. II 19961). For the most promptly exp 



1995 
oding 

systems, the companions may be 5-6 Mq main sequence 
(MS) sub-giants, with radii R± ~ 5 x 10 11 cm. More 
commonly, the MS companions may be 1 — 3 Mq sub- 
giants with radii ~ 1 - 3 x 10 11 cm. In the red giant 
(RG) case, the stars are evolved ^1 — 2 Mq stars with 
R* ~ 10 13 cm. In most scenarios, the companion is be- 
lieved to be in Roche lobe overflow. The separation dis- 
tance, a, is then comparable to the companion radius; 
for typical mass ratios, <z/iZ* = 2 — 3. 

After the supernova explodes, the ejected debris ex- 
pands freely for some time before hitting the companion. 
The flow becomes homologous, and the radius of a fluid 
element r — vt where v is the velocity and t is the time 
since explosion. We will d escribe the ejecta density pro- 
file by a broken power law ( Chevalier fc Soker|1989[ ) with 
a shallow inner region p\ oc r~ d and a steep outer region 



Po ex r n . The profiles join at the transition velocity 
vt = 6 X 10 s UE 51 /M C )^ 2 



cm s 



(1) 

where E^i = E/10 51 ergs is the explosion energy, M c — 
M/Mch is the ejecta mass in units of the Chandrasekhar 
mass, and is a numerical constant. The density in the 
outer layers (y > v t ) is 



Po(r,t) = C 



p v?t 3 \v t t 



r 



(2) 



with a similar expression for the inner layers. The nu- 
merical constants follow from the requirement that the 
density integrate to the specified mass and kinetic energy 



2(5-<T)(n-5)' 
(3-<J)(n-3) . 
1 (n-3)(3-£) 
Att n — 5 



1/2 



(3) 



The broken power law profile was originally derived for 
core collapse supernovae, but it fits multi-dimensional 
delayed-detonation models of SNe la remarkably well. 
For the models of Kasen et al. (20091 we find typical 
values of S = 1, n = 10, in which case the constants are 
& = 1.69, C P = 0.12. 

The characteristic timescale for the supernova ejecta 
to interact with the companion is 



U = a I vt 



10 4 ai3V 9 1 sec, 



(4) 

For 
5 hours, 



where ais = a/10 cm, and vg — ft/10 
RG companions the interaction timescale ti 
while for MS sub-giants, ss 5 — 30 minutes. 

The ejecta is highly supersonic when it collides with the 
companion, with mach number Ai ~ (a/i? w( j) 1 / 2 3> 1. 
Figure 1 illustrates the hydrodynamics of the interac- 
tion in a 2-di mensional numerica l calculation using the 
FLASH code ( Fryxell et al. 2000 1 and assuming a poly- 
tropic 7 = 4/3 equation of state, appropriate for a radi- 
ation dominated gas. As the flow sweeps over the com- 
panion star, a bow shock forms. Ejecta passing through 
the shock is heated and compressed into a thin shell, and 
its velocity vector is redirected. 

We will rely here on an approximate analytic descrip- 
tion of the collision dynamics. Material moving at ve- 
locity v interacts at a time t v ~ a/ v. The ejecta 
properties immediately after being shocked are given by 
the Rankine-Hugoniot jump conditions in the hypersonic 
limit. The density of the shocked gas is 



Ps{v) 



7 + 1 



7- 



1 



po(a,t v ) = 7£ 



M 

P ~~ T 



-n+3 



(5) 



taking 7 = 4/3. The pressure of the shocked gas is of 
order the incoming ram pressure 



Ps(v) 



1+7 



p v 2 sin 2 x 



C P sin 2 x- 



Mvf 



Vt 



-n+5 



(6) 

where x is the angle of the oblique shock front relative 
to the flow direction. The actual value of x varies along 
the bow shock, but for simplicity we take a constant, 
characteristic value near the maximal turning angle for 
hypersonic flows, x ~ 45°. 
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For a radiation dominated gas, p s = a^T 4 /3, where 
an is the radiation constant, which gives the equilibrium 
temperature of the shocked debris 



T s (v) = 2.8 x 10 6 MVVV3 3/4 (y 



-(n-5)/4 



K (7) 



For a RG companion at a = 2.5 x 10 13 cm, the outer 
layers of ejecta (v ~ 3vt) have T s pa 3 x 10 5 K, a value 
confirmed by the numerical calculation (Figure la). MS 
companions will have higher shock temperatures with 
T s i=s 10 6 — 10 7 K. To check the assumption of radiation 
domination in the shocked region, we note that Eqs. [5] 
and [6] imply a ratio of radiation energy to electron/ion 
energy 



aT 4 



, = 707 a^M, 

p s k B T s /m p 



3 / 4 iU-V4 3/2 



(8) 



which is ^> 1 for the scenarios under consideration. The 
details of thermalization, however, deserve further in- 
vestigation. Initially, the electrons/ions are shocked to 
very high temperatures (~ 10 9 — 10 10 K), then radia- 
tively cool toward the equilibrium value Eq. [7] by pro- 
cesses similar to those in ordinary supernova shocks - 
i.e., bremsstrahlung followed by Com pton up-scatte ring, 
and in some cases pair-production (Weaver 1976). If 
the timescale for the gas to radiate lags the dynami- 
cal timescale, non-equilibrium temperatures significantly 
greater than T s can be realized in the relaxation region. 
This is likely to occur in the highest-velocity, lowest- 
density ou termost layers, an d may be a source of harder 
radiation ( Katz et al.| 2009 1. 

The interaction with the companion diverts the incom- 
ing flow, carving out a conical hole in the supernova 
ejecta. The half opening angle of the hole is roughly 
#h = tan (rb/o) where r\> is the extant of the bow 
shock. Simulatio ns find r b ~ 2R+ and 9^ = 30 — 40° 
(see Figure 1 and |Marietta et ah (2000 1). The solid an- 
gle, Qh, of the hole is 



Oh 

47T 



cos 9h 



1 



1 



1 + (r b /a,y 



— ( 9 ) 
10 v ; 



This hole will provide a channel for radiation to quickly 
escape from the otherwise optically thick ejecta. 

The ejecta displaced from the hole piles up into a com- 
pressed shell along the cone surface. Assuming this shell 
layer is thin, its thickness Z s h can be estimated by mass 
conservation. The volume of a region of radial extent dr 
within the conical cavity is V\ = Q,\ l a 2 dr. The gas swept 
out of this region occupies a volume Vf = 2iral s ] 1 dr. The 
condition poVi — p s Vf gives 



l s h _ fih 2po ^ 1 
a 4ir p s 35 



(10) 



A dense shell of roughly this thickness is seen in Fig- 
ure la. The actual dynamics can become quite complex, 
with the shell broken in pieces by shear instabilities and 
the companion envelope shredded. 

After passing by the companion star, the shocked gas 
can expand laterally to try to refill the evacuated hole. 
The situation rese mbles the isentropic expansi on of a gas 
cloud into vacuum ( [Zel'Dovich fc Raizer|19 67 ) ; the front 
of the rarefaction wave moves outward at the maximum 



escape velocity vi = 2/(j—l)c s where c s = {^pg/ ps) 1 / 2 — 
(8/49) 1 / 2 sin x v is the sound speed of the shocked ma- 
terial. The net velocity in the direction perpendicular 
to the symmetry axis is then v x — v sin 9h — vi cos 9^ rs 
—0.2 v. The ejecta moving at velocity v thus re-closes at 
a time Fi k /\v x \ « 5i?*/u after passing by the companion, 
or at a time th = (a + 6R+)/v after the explosion. In the 
adiabatic calculation (Figure lb) , the rarefaction softens 
the density gradient in the polar direction, but fails to re- 
fill the shadowcone region uniformly before freezing out. 
If radiative cooling is significant during these phases (see 
Section [3]) the sound speed will be reduced, which would 
further delay or halt the closing of the hole. 

The energy density of the shocked gas is e s = 3p s and 
so the total energy dissipated in the collision shock is 
found to be 



. is o „ 

^49 Sm *4^ 



1.5 x 10 49 E 51 ergs. 



(11) 



Much of this thermal energy will be lost again to adi- 
abatic expansion, but if even a fraction is radiated the 
collision luminosity should be quite bright. 

3. PROMPT X-RAY BURST 

As supernova ejecta flows past the companion star, the 
hot surface layers of the shocked shell become exposed 
(Figure la). At this time, some radiation may be able to 
escape straightaway through the evacuated shadowcone 
hole, giving rise to a sudden burst of emission. In general, 
only a fraction of the total energy in the shell can be 
radiated promptly - i.e., before suffering significant loses 
due to adiabatic expansion. This prompt emission arises 
from the surface layer of the shell with a thickness, la, 
determined by requiring that the diffusion time through 
that layer 



td = T— = 

3c 



3c 



(12) 



be less than the timescale for expansion, given by 
the shell sound crossing time l s h/c s (which is typically 
shorter than the dynamical timescale a/v). Using Eqs.p] 
and [lO] gives 



Zd a ( 4-7T 
where the quantity 



1/2 



14C P sin x 



1/2 



(n-4)/2 



( kM\ 



1/2 



29 My V /2 «e /2 days 



(13) 



(14) 



1982) 



is the familiar "effective" diffusion time (Arnett 
that sets the duration of the ordinary 56 Ni-powered SIN la 
light curve. We have assumed a constant opacity K e = 
0.2 cm 2 g _1 , appropriate for electron scattering in fully 
ionized A/Z = 2 elements. For RG companions at a w 
10 13 cm, we find id ~ Z s h for the layers v > 2u t . In this 
case, most of the energy dissipated in the outers layers 
can be radiated in the prompt burst. For MS companions 
with a ~ 10 11 - 10 12 cm, the ratio l d /l sh ~ 0.01 - 0.1, 
and only a fraction of the photons escape promptly. 

While the bulk of the SN ejecta remains extremely op- 
tically thick at this phase, photons can initially escape 
through the hole carved out in the interaction. This 
channel will close, however, once the outer layers of ejecta 
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have re-expanded to fill the shadowcone, which happens 
at a time th ~ (a + 6i?*)/v max where w max is the max- 
imum ejecta velocity (see Section [2|. A given layer of 
ejecta can contribute to the burst only if it passes the 
companion at a time less than th, which holds for mate- 
rial moving faster than -y m ; n ~ w max (a + R+)/(a + 6i?*). 
For a/R* = 3, we find ?; m i n « u max /2. 

By integrating the energy density e s = 3p s , of the 
shocked gas within l d , we can estimate the total energy 
escaping in the burst 



£ x = 2tt 



3p s (v)l d (v) a - )dv 



7.72 x 10 4 



-in 



C P sin 13 x 



1/2 



G2 



(15) 



013 



(n-6)/2 



ergs, 



where we approximated the upper limit as K max — > oo. 
The duration of the burst is the time it takes the 
emitting ejecta to flow past the companion, or At x — 
(a + -R*)/Vmi n — (a + R+)/v mei x- For typical values 
(v m in ~ 2«t; u ma x ~ 4t> t ; i?* » 3a) this duration is 
Ai x ti/3. Assuming the radiation is emitted into 
a solid angle f2h, the isotropic equivalent luminosity is 
L x iso = (47r/f2h)(-Ex/Af x ). Taking characteristic values 
(n,'i,xA,»min) = (10,l,45°,40°,2 Wt ) we find 



.8 x 10 4 



ergs s 



(16) 



This luminosity is independent of a, and so roughly com- 
parable for all types of companions. The value is similar 
to that of shock breakout in core collapse SNe, which 
is not surprising given that the shock temperature and 
emitting surface area are comparable in the two phenom- 
ena. The collision burst will only be visible for viewing 
angles peering down the hole, 6 < 9^. Such an orienta- 
tion occurs f2h/47r as 10% of the time. 

The spectrum of the prompt burst may be approxi- 
mated as a blackbody at the equilibrium shock tempera- 
ture (Eq. [7| , implying emission peaking in the soft X-ray 
with typical energies T x ~ 0.05 — 0.1 keV for RG and 
T x ~ 0.2 — 2 keV for MS companions. Non-equilibrium 
effects (see Section |2| could lead to some emission at 
significantly higher energies (10 — 100 keV), while non- 
thermal particle acceleration may also contribute a power 
law continuum of hard radiation. X-rays emitted in the 
direction of the companion star will ionize its surface 
layers and be reprocessed, likely giving rise to sub stan- 
tial line recombination/fluores cence emission (e.g., Bal- 
lantyne & Ramirez-Ruiz 20011) 



According to Fq. [7] the temperature of the burst spec- 
trum evolves with time as 



r„(t) = o.i «r 3 3/ X 1/ % 1/2 (^) 



(n-5)/4 



keV, (17) 



which predicts that, at least initially, the spectrum be- 
comes harder with time, as interior layers of ejecta have 
higher densities and shock temperatures. On the other 
hand, deviations from equilibrium arc likely to be great- 
est in the highest velocity layers, which may counter- 
act this trend. Eventually, as ejecta gradually refills the 
shadowcone, the effective photosphere moves to a larger 
radius and the emission must soften. 
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Fig. 2. — Model light curve of a Type la supernova having 
collided with a red giant companion at a separation distance a = 
2 X 10 13 cm. The luminosity due to the collision is prominent 
at times t < 8 days and for viewing angles looking down on the 
collision region (8 = 0°). At later times the emission is powered 
by the radioactive decay of 0.6 Mq of 56 Ni located in the inner 
layers of ejecta (v < 10 9 cm s" 1 ). The black dashed line s how s 
the analytic prediction for the early phase luminosity (Eq. |22| . 

The actual structure of the collision region will be more 
complex and inhomogeneous than that imagined here, 
due either to the inherent dumpiness of the supernova 
ejecta, or to secondary shocks and hydrod ynamical insta- 
bilities dev eloping in the interaction (e.g 



et al. 1996 



Cid-Fernandes 
This may lead to fluctuations in the burst 
light curve on a time scale SR/v where SR is the typical 
clump size. Multi-dimensional radiation- hydrodynamics 
calculations will be needed to characterize the light curve 
and spectra in detail. 

4. EARLY UV/OPTICAL LUMINOSITY 

Thermal energy not radiated in the prompt burst can 
diffuse out in the hours and days that follow, but will 
suffer loses from adiabatic expansion. At t ~ 1 day, this 
emission will be primarily at UV/optical wavelengths. 

We consider here times ^> t\ such that homology has 
been re-established in the ejecta. The final ejecta struc- 
ture will be asymmetric, but for now we neglect an- 
gular dependencies. The density profile is taken to be 
p{{r,t) — 7/o po(r, i), where the constant /q < 1 ac- 
counts for the lateral expansion and radial readjustment 
that occur during the transition to homology. 

Assuming the evolution is adiabatic, the pressure pro- 
file of a shocked fluid element at these times is pt(t) = 
p s [Pf(t) / PsV ■ For the region affected by the collision 
(6 < 0h) 



6 



Vt{r,t) = y/o /3 CpSin 2 x 



Ma 



v(t A \ v t t J 



-n+l 



(18) 



The pressure is negligible outside 9^. As time progresses, 
the adiabatic profile Eq. [18] will continue to describe the 
opaque inner layers of ejecta, but the outer layers will 



Supernova Companion Collision Emission 



5 



be modified by radiative diffusion. The evolution can 
be calculated us ing a self-similar diffusion wave analysis 
( Chevalier|f992 1. At a time t, diffusion will have affected 
the ejecta above a radius r<i determined by setting the 
diffusion time td — r^npf/3c equal to the elapsed time t 



Cp ,2 

3 sn 



1 l/("-2) 



t (n-4)/(n-2) 



(19) 



Over time, the diffusion wave recedes into the ejecta in 
a Lagrangian sense. However, for times f < 5 days, 
remains in the steep outer layers of ejecta (v > i>t). 

Considering only the radial transport, the isotropic 
equivalent luminosity in the comoving frame is given by 
the diffusion approximation 



-^c.iso ^) 



apt dr 



(20) 



We continue to assume a constant opacity, even though 
at UV wavelengths the line expansion opacity may exceed 
electron scattering and will, to some extant, be temper- 
ature dependent. 

In the outer layers of ejecta (r > r^) the comoving 
luminosity L Ci ; so is constant with radius. Its value is 
therefore set by processes near the diffusion wave radius. 
A reasonable estimate of -L Cj i so is derived by evaluating 



Eq. 20 at ra, using the pressure profile Eq. 18 



L c ,iso = ~ l)?7sin 2 x/d /3 



2/(n-2) 



-4/(n-2) 



(21) 



where 77 is a constant of order unity that must be deter- 
mined by solving the full diffusion equation (Chevalier 
19921. Written in this form, the resemblance to core 
collapse supernova light curves is clear: apart from con- 
stants, the luminosity is of order E/t sn times a factor 
ti/tsn that accounts for adiabatic loses. 

Taking (n, 6, x,fo,v) = (10, 1, 45°, 0.5, 0.5) we find that 
from appropriate viewing angles 



10 43 013 M, 



1/4,7/4-3/4,-1/2 



W er g ss " 



(22) 



where tday is the time since explosion measured in days. 
The derivation applies only at times significantly greater 
than ij, but Eq. [22] may provide a workable estimate at 
earlier times. The observer frame luminosity differs from 
Eq.[22]by an additional term accounting for the advected 
luminosity, of order v/c < 10% compared to L c ,iso- 

The collision luminosity will only be readily discernible 
when it exceeds the luminosity, L n i, of the ordinary 
56 Ni powered light curve. At early times, t <C t sn , 
the approximate analytic light curves of SNe la give 
L ni = e ni M ni (t/t sn ) 2 , where e ni = 4.8 x 10 10 ergs s^g" 1 
and M ni is the mass of 56 Ni ( |Arnett|l982| . We then find 



that L r 



> L n [ for times 



t < 7.3 a% b M, 



2/5 



(23) 

where M n ; .0.6 = M n ;/0.6M Q . Note that the opacity in 
the 56 Ni region, K n j, is heavily affected by iron group 
line blanketing, and so may be greater than the opacity 



n e in the outer layers. This would help delay the 56 Ni 
luminosity takeover. 

The wavelength of the early emission depends on the 
photospheric radius, r p , defined as the location where the 
optical depth r = ptndr = 1, or 



f (n-3)/(n-l) ( CpKMv, 



(n-3) x l/(n-l) 



1 

The effective temperature of the emission 
(^c,iso/47rrpCr) 1 / 4 is then, for n = 10 

T„ff = 2.5 x 10 4 V 35/36 t 



(24) 



-37/72 



K. 



May - (25) 

At t = 1 day, the emission peaks at wavelengths A ss 
1000 A, however the Rayleigh- Jeans tail of the blackbody 
extends into the near UV and optical. 

While the analytic solution captures the essential phys- 
ical effects, it neglects the ejecta asymmetry and non- 
radial transport. We have therefore calculated numer- 
ical light curves using a 3-D non-grey radiation trans- 
fer code which also ac counts for the effec ts of line opac- 
ity and 56 Ni heating ( |Kasen et al. 2006| . For simplic- 
ity, we use an artificial ejecta model with density profile 
pt(r,0,t) — p (r,t)fo(8), where the function fo(0) now 
describes the angular structure of the conical hole region 

h(0) = fh + (1 - fh)r^. [l + Aexp 



1 



(o P /o h y 



( 26 ) 

where x = 9/6^. This formula approximates the results 
of simulation in the homologous phase choosing m = 8, 
f h = 0.1. 6» h = 30°, 9 p = 15°, and A = 1.8. The pressure 
profile in the shocked region was taken from Eq. 18] 

The light curve calculation (Figure 2) shows that the 
early collision luminosity is dramatic when the compan- 
ion is a RG at a — 2 x 10 13 cm. The numerical results 
agree with the analytic estimate (Eq. [22| , and also illus- 
trate the anisotropy of the radiation. The collision lumi- 
nosity is brightest for viewing angles looking down upon 
the shocked region (6 < 8^), but a significant amount 
of radiation diffuses out at angles 9 w 90°, and a few 
percent is even back-scattered along 9 sa 180°. 

While the companion interaction produces a conspicu- 
ous signature (a kink) in the early bolometric light curve, 
most of the emission is in the UV; in the optical bands, 
the effect is less dramatic (Figure 3). For a RG com- 
panion, the B-band light curve shows a distinct bump 
at t < 5 days which should be relatively easy to probe 
observationally. At redder wavelengths, or for smaller 
separation distance a, the collision simply modifies the 
shape of the light curve rise. However, because SN la 
light curves are quite standard, a statistical analysis of 
(good quality) early time photometry should be able to 
pull out these subtle differences. 

5. OBSERVATIONAL PROSPECTS 

The results derived here suggest a new means for 
constraining supernova companions using early photo- 
metric observations. Table 1 summarizes the analytic 
estimates of the collision emission for various SNe la 
progenitors. The theoretically predicted signatures ap- 
pear to be quite robust, as they rely only on estab- 
lished physics familiar from the core collapse SNe con- 
text. However, further numerical studies using multi- 
dimensional radiation-hydrodynamical calculations (and 
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TABLE 1 

Properties of Type Ia Supernova Collision Emission 



Analytic Estimates 



Companion 


a 








At x 


-^x,iso 


T x 


^ciso (1 day) 


tc 


RG (M ~ 1 M ) 


2 x 10 13 


cm 


3.9 x 10 47 


ergs 


1.9 hours 


5.8 x 10 44 


0.07 keV 


2 X 10 43 


9.6 days 


MS (M ~ 6 Mq) 


2 x 10 12 


cm 


3.9 x 10 46 


ergs 


11 mins 


5.8 x 10 44 


0.2 keV 


2 X 10 42 


3.8 days 


MS (M ~ 2 M ) 


5 x 10 11 


cm 


9.6 x 10 45 


ergs 


2.8 mins 


5.8 x 10 44 


1.0 keV 


5 X 10 41 


2.2 days 


MS (M ~ 1 Mq) 


3 x 10 11 


cm 


5.8 x 10 45 


ergs 


1.7 mins 


5.8 x 10 44 


1.4 keV 


3 x 10 41 


1.8 days 





days since explosion 



days since explosion 



Fig. 3. — Signatures of companion interaction in the early 
broadband light curves of Type Ia supernovae. We model three 
different progenitor scenarios: a RG companion at a = 2 X 10 13 cm 
(green lines), a 6 Mq MS companion at a = 2 X 10 12 cm (blue 
lines) and a 2 Mq MS companion at a = 5 X 10 11 cm (red lines). 
The ultraviolet light curves are constructed by integrating the flux 
in the region 1000 — 3000 A and converting to the AB magnitude 
system. For all light curves shown, the viewing angle is 6 = 0° . 

including non-equilibrium effects) will be needed to refine 
the detailed light curve and spectrum predictions. 

For all companion types, signatures of the collision will 
be prominent only for viewing angles looking down upon 
the shocked region, or ~ 10% of the time. Detection 
will therefore require high cadence observations of many 
supernovae at the earliest phases (< 5 days) and at the 
bluest wavelengths possible. Ironically, these observa- 
tions may sometimes be easier for distant SNe. At red- 
shifts z > 0.5, the UV flux would be redshifted into the 
U-band, while cosmological time dilation would prolong 
the light curve by a factor (1 + z). 

Detecting the collision signatures becomes significantly 
easier for larger separation distances. Current optical 
and UV data sets likely already constrain red giant com- 
panions (a ~ 10 13 cm). Ongoing or upcoming surveys 
could be tuned to probe the larger (M > 3 Mq) main 
sequence companions (a ~ 10 12 cm). Optical detec- 
tion of the smallest ~ 1 Mq main sequence companions 
(a ~ 10 11 cm) will be challenging, requiring measure- 
ment of subtle differences in the light curves at t < 2 day. 
However in all cases the prompt X-ray burst sh ould be 
bright. P roposed X-ray surveys (e.g., EXIST, |Grind 



lay 20051 may then detect a large number of collision 
bursts every year, at least in the case of MS companions 
which produce harder radiation. If non-equilibrium or 
non-thermal shock effects are operative, some hard radi- 
ation may accompany all bursts. 



The most compelling reason for studying the collision 
emission is that it offers a straight-forward measure of 
the separation distance between the stars. This value 
can be determined from the duration of the X-ray burst 
(At x ~ a/v), or its temperature (Eq. [7|, or from the 
luminosity of the early optical/UV emission (Eq. 22 1. If 
we in turn assume that the companion is in Roche lobe 
overflow, it's radius can be inferred i?* ~ a/3 — a/2. 
In principle, the ratio a/i?* itself could be constrained 
using a statistical sample of SNe Ia, as the anisotropy 
of the emission depends on the opening angle, #h, of the 
shocked region of eject a. 

While our discussion has focused on Type Ia super- 
novae, similar signatures of companion interaction should 
apply to Type Ib/Ic and some Ty pe II SNe which may 
arise from close massive binaries ( |Podsiadlowski et al. 
19921. In these cases, the shock breakout from the ex- 



ploding star will contribute to the luminosity on compa- 
rable timescales, likely producing two peaks in the X-ray 
emission. Some gamma-ray bursts might also come from 
binary systems, and the interaction of the relativistic ma- 
terial with a st ellar companion may pro duce another type 
of X-ray flare ( |MacFadyen et al.||2005| ). 

It is also possible that early emission from SNe stems 
from a collision not with the companion star, but with a 
surrounding circumstellar medium (CSM). To substan- 
tially decelerate the ejecta, the CSM would need to have 



113 , 



amass ~ 0.01-0.1 M Q located at radii ~ 10 11 — 10 1 
A slow (10 km s _1 ) stellar wind moves beyond these dis- 
tances is less than a year, so it may be difficult to realize 
these conditions in a single degenerate scenario of SNe Ia. 
In the double degenerate merger scenario, the total mass 
of the system can exceed Af c h, and a few 0.1 Mq of excess 
carbon/oxygen may linger in the vicinity. If this material 
remains at the tidal radius ~ 10 9 cm, the resulting emis- 
sion will be extremely brief (~ 1 sec), however if some 
mass is puffed out to larger radii in the super-Eddington 
accretion phase of the merger, the emission may be sim- 
ilar to that discussed here. Interaction with a spherical 
CSM is distinguishable from companion interaction by 
its luminosity function; in the former case, the emission 
should be nearly the same from all viewing angles. 

In either case, the early time emission of supernovae 
provides much needed insight into the nature of the pro- 
genitor system. Observational surveys could be designed, 
either from space or the ground, to acquire the collision 
signatures in a systematic way. If one collects a signif- 
icant number of events, it will be possible to correlate 
the measured separation distances with the properties of 
the ordinary 56 Ni powered light curve and spectra. Such 
observations would provide direct, empirical insight into 
how the parameters of the progenitor system influence 
the outcome of supernova explosions. 
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